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Abstract 

In the development of the first generation of sensors and memory chips 
based on spin-dependent tunneling through a thin trilayer, it has become 
clear that pinhole defects can have a deleterious effect on magnetoresistance. 
However, current diagnostic protocols based on Andreev reflection and the 
temperature dependence of junction resistance may not be suitable for pro- 
duction quality control. We show that the current density in a tunnel junc- 
tion in the cross-strip geometry becomes very inhomogeneous in the presence 
of a single pinhole, yielding a four-terminal resistance that depends on the 
location of the pinhole in the junction. Taking advantage of this position de- 
pendence, we propose a simple protocol of four four-terminal measurements. 
Solving an inverse problem, we can diagnose the presence of a pinhole and 
estimate its position and resistance. 

PACS '03: 73.40.Rw, 73.40. Jn, 85.75.Dd, 85. 75. Mm 



Semiconductor manufacturers are currently developing magnetic-random-access-memory 
elements 1-3 based on magnetic tunnel junctions; such junctions separate two ferromagnetic 
metallic leads by a thin insulting layer, 4-6 often made by oxidizing a film of Al or other 
suitable metal. Both because of the thinness of the insulating layer and because of the 
possibility of inadequate oxidation, "pinhole" defects — direct metal-metal shorts through 
the nominal insulator — have attracted significant attention. 7-13 A single pinhole can also 
be generated in a previously pinhole-free junction carefully through a voltage ramp 14-17 
or, by implication, inadvertently. Generally, the parasitic current through pinholes detracts 
from a junction's magnetoresistance, 18 so methods for diagnosing and locating such defects 
become important during the development of practical devices. Surprisingly, a fit of differ- 
ential conductance to the Simmons form 19 fails unambiguously to guarantee the absence of 
pinholes. 11,20 Surer methods include the use of an integrated superconducting electrode, 20 
the temperature dependence of device resistance, 20-22 and surface decoration. 7 ' 14 

These diagnostics may not be integrated easily into a development or manufacturing 
process. We therefore propose a very simple test that not only determines the presence 
of a pinhole with high confidence but also can typically locate the pinhole to within 7% 
of the junction area and estimate the pinhole resistance. We propose four four-terminal 
measurements. A discrete three-dimensional resistor model allows us to compute the result 
of each measurement for an assumed pinhole position; working backward from a set of four 
experimental measurements to the pinhole position therefore constitutes an inverse problem, 
to which we demonstrate a solution. Since the method applies to non-magnetic as well as 
to magnetic tunnel junctions, we shall generally treat the junctions without reference to 
magnetic properties. 



We will consider the common cross-strip geometry of figure 1, consisting of a conduct- 
ing top strip separated from a perpendicular conducting bottom strip by a thin insulator. 
The device's four leads are numbered as in the figure. In a typical four-terminal resis- 
tance measurement, as shown, current is injected in lead 1 (top strip) and removed through 
lead 4 (bottom strip), while the voltage is measured between leads 3 and 2. Moodera and 
collaborators 23-24 have pointed out that, because a conducting strip is not an isopotential, 
such a four-terminal measurement will give misleading, sometimes even negative, absolute 
resistances when the resistance of the insulating layer is comparable to or smaller than those 
of the upper and lower metals. We generalize their observation for a junction shorted by a 
single pinhole. 
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Figure 1. The standard cross-strip geometry (exploded view) consists of a lower metallic 
strip over which is deposited an aluminum layer, which is oxidized (shaded layer) before 
a top metallic strip is deposited. To obtain the four-terminal resistance R\, one injects 
current through the leads labeled 1 and 4 while measuring the voltage between leads 3 and 
2. A pinhole short through the insulating layer will result in an unreliable four-terminal 
resistance that depends on the position of the pinhole. 

We model the junction as a discretized three-layer resistor network. Each of the four 
long leads is assumed isopotential where it makes contact with the square junction. (Later, 
we address this assumption as a source of uncertainty.) The pinhole is modeled as a metallic 
(Al) inclusion of diameter approximately 0.8 nm penetrating the insulating layer. Classical 
conductivities are used for conduction through metallic channels, including the pinhole, while 
the very low conductance of tunneling is approximated from Simmons's formula assuming a 
barrier height about 2 eV and a thickness around 1 nm. Kirchhoff 's law yields a set of linear 
equations for the potentials at all the nodes of the resistor network in terms of fixed current 
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/ through two leads. We solve the linear equations numerically in order to calculate the 
voltage difference V between the two measurement leads. 

Denote pinhole resistance by R p and the intrinsic tunnel "resistance," determined in the 
Simmons model by the barrier height, area, and thickness, by Rt. Figure 2 plots the nominal 
four-terminal resistance (R\) against Rt for R p fixed at 50 Q (horizontal line); the top and 
bottom layers are given two-dimensional resistivities of R io „ = 10Q/D and R hot = 20Q/D. 
These values appear representative of recent experiments. 2 ^-24 The pinhole for this figure 
was placed in the bottom-left quadrant but not very near the corner. As noted, the nominal 
four-terminal resistance takes anomalously small, even negative, values when Rt -C -R top , -Rbot- 
The four-terminal measurement is independent of Rt when the latter is large compared to 
the resistances of the leads, since then most current flows through the pinhole. One thus 
measures a four-terminal resistance 

R 1 (r) = R p + R lf (r) , (1) 

where R\f is a function of the pinhole's position, r, for fixed R top and R bot . In a junction of 
sufficiently small area, i? top , R bot , Rp <C Rt, leaving the positional term, Rij(r), independent 
of R p . The inset to figure 2 shows voltage contours on the top layer; the pinhole is evident. 
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Figure 2. The four-terminal resistance of a junction with a pinhole is plotted as a function 
of the effective resistance of the tunneling layer. For this illustration, the pinhole resistance 
Rp = 50f2, while the top- and bottom-layer resistivities take the values 10 O/D and 20 fi/D 
When the tunneling resistance is sufficiently large, most of the current flows through the 
pinhole, and the four-terminal resistance is the sum of the pinhole resistance R p and a 
position-dependent piece i?ij(r), as in (1). The inset shows constant-voltage contours on 
the top layer. 

We now iterate over all possible pinhole positions, r. Figure 3 illustrates, for the same 
parameters as figure 2, the position dependence of -Ri(r). If we knew the value of R p , this 
figure would localize the pinhole to the curve of constant R\ corresponding to the measured 
value. 
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Figure 3. The predicted four-terminal resistance R\ is plotted against pinhole position. 
By rotating the roles played by the four leads, one can localize the pinhole over most of 
the junction area. 

In fact, we can find both R p and r simply by rotating the roles played by the four 
electrodes. Thus we measure R 2 by forcing current through leads 1 and 2, measuring voltage 
across leads 3 and 4. Repeating the rotation, we get four four-terminal measurements, as in 
table I. In the absence of a pinhole, these four measurements would be equal. 





I 


V 


Ri 


1,4 


3,2 


R2 


2, 1 


4,3 




3,2 


1,4 


R A 


4,3 


2, 1 



Table I. Four four-terminal resistances (R\ ■ ■ ■ R4) are defined from the four ways of 
injecting current (I) through one top and one bottom lead while measuring the voltage 
(V) between the remaining two leads. 



From (1), differences of the resistances in table I will not depend on the pinhole resistance 
R p ; it is convenient 25 to take as a maximal independent set 

R a = R 1 -R 2 + Rs-Ra 

R b = Ri + R 2 -Rz-Ra (2) 

R c = R l -R 2 -R 3 + R A . 
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Figure 4 displays contour maps of each of the resistances from (2) against position r of 
the pinhole, with the contour corresponding to the (simulated) measurement outlined. If we 
may assume perfect knowledge of the geometry, the three outlined contours will intersect at 
a point, thus locating the pinhole. 

Measurements on an actual tunnel junction will suffer errors and uncertainties from 
several sources, such as surface roughness and possible non-percolating conducting inclusions. 
The classical model itself breaks down for junctions small compared to the relevant mean 
free path. We model the uncertainties in two ways, chosen to be representative rather 
than microscopically realistic. We then propagate the uncertainties to the calculation of 
resistances as a function of pinhole position. 

First, we consider a break of up to 10% of the junction width at the point where one lead 
joins the junction, as though there were a bad solder joint. Of course, there is no solder — 
the strip is continuous — but we use this 10% void as a proxy for uncertainty in the actual 
junction geometry. A full calculation would find equipotentials bending into the leads where 
the current is injected rather than coming in straight; inhomogeneity in the leads would 
further complicate the picture. The 10% void is varied over positions at this junction edge, 
leading to an error in the four-terminal measurement (compared to the result without the 
break) . 

Second, we posit a single atomic-scale terrace in the insulating layer; this too is allowed 
to migrate over the junction, leading to an error. This source of error stands in as well for 
possible tunneling "hot spots" due to roughness, 26-29 so long as the tunnel current remains 
small compared to the pinhole current. The two sources of uncertainty are combined to give 
estimated error bars in the resistances R a , R^, and R c . These error bars can be pictured as 
broadening the outlined contours of figure 4, and the single point of overlap becomes instead 
a region in which we expect to find the pinhole. 




Figure 4. Resistance contours for R a , Rf,, and R c from equation (2). The resistance that 
would be measured for the example of figure 2 is emphasized in each. The pinhole is located 
at the intersection of the three emphasized contours. Because of the two (perpendicular) 
node lines in R a , the values are much smaller: the contour scale on R a is smaller by a 
factor of 20 than the scales for the other two resistances. 



For the sample geometry we have been considering, these sources of error yield uncer- 
tainties in each of the four-terminal measurements of 0.5%. If all four measurements Ri 
(equation (1)) agree within this percentage, 30 we conclude with high confidence that the 
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junction does not harbor a pinhole. (We discuss below the probability of a false negative 
result.) 

We have discussed the calculation of R- t given a known pinhole, but in fact we wish to 
determine the position of a pinhole given only the measurements R{. We solve this inverse 
problem by brute-force computation of four-terminal resistances for all possible positions of 
a single pinhole on a two-dimensional grid. Confidence regions, as in figure 5, are determined 
by x 2 minimization and verified with Monte-Carlo sampling. 31 
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Figure 5. The pinhole of the example localized. The two contours bound 68% and 95% 
confidence regions as calculated with the x 2 method. (Monte-Carlo gives similar pictures.) 

To extract a confidence region from the x 2 method, we compare experimental measure- 
ments to calculated resistances Ri(r) as functions of pinhole position. The objective 
function x is the sum of square errors between calculated and measures resistances, nor- 
malized by the estimated variances a 2 : 



X 2 (r) 



E (hw-^ (0) ]M) 



(3) 



-a.b.c 



We apply the usual chi-square distribution with two degrees of freedom (corresponding to 
r) to convert x contours into confidence regions. 

To verify these results with Monte-Carlo sampling, we start with an assumed pinhole 

position and resistance and calculate the four resistances R.f^ ■ We then generate a large 

number of "experimental" data sets R^ by adding Gaussian-distributed noise according 
to the standard deviations worked out previously. For each data set, we find the pinhole 

position r (on the discrete grid) that minimizes x 2 - With a large number of sets R^\ we 
accumulate the number of times each grid position minimizes the objective function. In the 
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68% confidence region, we include first the grid point with the largest count, then that with 
the second- largest count etc. until we have included 68% of the data sets. The results agree 
closely with those from the x analysis. The procedure could be refined by increasing the 
number of grid points in the vicinity of the most probable pinhole position until the resolved 
size fell below the dimensions of the 68% (1-sigma) confidence bound. 

We have repeated this calculation for pinhole positions covering the junction area. As 
the four- fold symmetry makes evident, a pinhole in the very center of the junction would be 
undetectable through this method, since all four (equation (1)) would be equal. Allowing 
for estimated errors in the four measurements, we find a region covering about 13% of the area 
around the center in which a pinhole cannot be distinguished from the absence of a pinhole: 
a pinhole in this central region still yields measurements R a , R^, R c all equal to zero within 
an uncertainty 2a. Confidence regions for pinholes near the central region are larger than 
for those closer to the edges and corners. Considering all possible actual pinhole positions 
(including the center), we find that, on average, the 68% confidence region comprises 3.5% of 
the junction area, while the 95% confidence region comprises 7%. Having located the pinhole, 
we can solve for R p in (1) or, to take advantage of averaging, in R^ = R\ + i?2 + R3 + Ra- 

We have proposed a simple technique for diagnosing, with high confidence, the pres- 
ence or absence of a single pinhole in a tunnel junction and furthermore for determining its 
position. We hope soon for an experimental test. Pinhole- free junctions may be prepared 
(and verified using the temperature dependence of resistance 22 ' 20 or absence of Andreev 
reflection 20 ) and subjected to the protocol outlined here. The four four-terminal measure- 
ments should be equal to within a small uncertainty. To the extent that this uncertainty, 
derived from a number of junctions, differs from our rough estimate, it will serve to recali- 
brate our estimated error bars. Then, applying a slow voltage ramp, a single pinhole can be 
created in each junction. The protocol should determine the pinhole's presence. If possible, 
it will be interesting to check actual pinhole location with a decoration method 7 ' 14 or by 
scanning-probe microscopy. 26,32 

Work now in preparation investigates the effects of Ohmic heating and thermal transport 
on the differential conductance of a junction incorporating both a pinhole and tunneling 
channels. Further work will replace the classical approximations made so far with quantum- 
mechanical calculations. 33 
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